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We show how pattern formation in Faraday waves may be manipulated by varying the harmonic 
content of the periodic forcing function. Our approach relies on the crucial influence of resonant 
triad interactions coupling pairs of critical standing wave modes with damped, spatio-temporally 
resonant modes. Under the assumption of weak damping and forcing, we perform a symmetry- 
based analysis that reveals the damped modes most relevant for pattern selection, and how the 
strength of the corresponding triad interactions depends on the forcing frequencies, amplitudes, 
and phases. In many cases, the further assumption of Hamiltonian structure in the inviscid limit 
determines whether the given triad interaction has an enhancing or suppressing effect on related 
patterns. Surprisingly, even for forcing functions with arbitrarily many frequency components, there 
are at most five frequencies that affect each of the important triad interactions at leading order. 
The relative phases of those forcing components play a key role, sometimes making the difference 
between an enhancing and suppressing effect. In numerical examples, we examine the validity of 
our results for larger values of the damping and forcing. Finally, we apply our findings to one- 
dimensional periodic patterns obtained with impulsive forcing and to two-dimensional superlattice 
patterns and quasipatterns obtained with multi-frequency forcing. 

PACS numbers: 05.45.-a, 47.35. +i, 47.54.+r, 89.75.Kd 
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I. INTRODUCTION 



Parametrically forced surface waves have proven to be 
a rich and versatile source of patterns since their ini- 
tial observation by Michael Faraday in 1831 [J. These 
Faraday wave patterns are composed of standing waves 
set up in response to periodic vertical vibration of suffi- 
cient strength. Early investigations (see 0j@l ^^^ reviews) 
used a sinusoidal forcing function and focused on simple 
patterns such as stripes, squares, and hexagons, which 
oscillate in subharmonic response to the forcing. Re- 
cently, experimentalists have used multi-frequency forc- 
ing to generate more complex states such as quasipat- 
terns and superlattice patterns 4_, _5,, ^ J^ ^, ^, J^ JT] . 
These observations have, in turn, fueled theoretical in- 
terest in such patterns and in multi-frequency forcing 
[H m El 15, 16, 17, ,ia, 19] . 

The use of multi-frequency forcing requires the se- 
lection of a large number of control parameters. The 
forcing frequencies, their amplitudes, and their relative 
phases may all affect the pattern formation problem in 
a nontrivial way jj, [12. ll4L |20). Further complexity 
arises from the presence of multiple length scales. In 
addition to the length scales driven by the various forc- 
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ing frequencies in accordance with the parametric {i.e., 
subharmonic) resonance conditions, there are numerous 
damped modes that are driven nonlinearly. Many of 
these can significantly influence the dynamics of the crit- 
ical modes. For instance, in the case of two-frequency 
forcing, the damped mode that oscillates at the so-called 
difference frequency is important for selecting superlat- 
tice patterns p^ . 

Resonant triad interactions - the lowest order nonlin- 
ear interactions - provide a useful framework for inves- 
tigating the relationship between the many control pa- 
rameters and length scales in the multi-frequency forced 
Faraday wave problem. Resonant triads that couple two 
critical modes with a damped, spatiotemporally resonant 
mode play a key role in the nonlinear pattern selection 
process. Most of these damped modes function as en- 
ergy sinks, effectively creating an anti-selection mecha- 
nism that suppresses the triad interaction and thereby fa- 
vors patterns which avoid the corresponding resonant an- 
gle. However, other damped modes act as energy sources, 
providing a positive selection mechanism that helps sta- 
bilize patterns involving the associated resonant angle. 
The effect of different damped modes on pattern selection 
is investigated in 20], which, for forcing functions with 
up to three frequency components, determines the most 
important damped modes, their effect (enhancing or sup- 
pressing) on associated patterns, and the dependence of 
the nonlinear interaction on the forcing frequencies, am- 
plitudes, and relative phases. These results are used to 



interpret recent Faraday wave experiments that produced 
complex patterns, namely, a two- frequency forced super- 
lattice pattern in [3| and a three-frequency forced quasi- 
pattern in [l3|. The approach developed in ^2U\ follows 
from a systematic consideration of weakly broken sym- 
metries: time translation, time reversal, and Hamilto- 
nian structure (see Il2l|), and is therefore most relevant 
for systems with weak damping and forcing. In this limit 
the vastness of the control parameter space can be an as- 
set, enabling one to enhance or suppress particular triad 
interactions simply by tuning the appropriate forcing pa- 
rameters. 

In this paper, we adopt the same prescriptive approach 
to Faraday wave pattern formation, describing in more 
detail the technique for exploiting weakly broken symme- 
tries, and extending the results of 0| to forcing functions 
containing arbitrarily many Fourier components. We de- 
termine which damped modes are favored by a strong 
nonlinear coupling and tabulate how the corresponding 
resonant triad interactions depend on the forcing param- 
eters. A somewhat surprising result, which makes this 
project feasible, is that for a given damped mode there 
are at most five out of the potentially infinite number 
of forcing frequency components in the forcing function 
that affect the resonant triad interaction at leading order 
in the damping parameter 7 (defined below). We inves- 
tigate numerically the validity of our predictions with 
respect to the small 7 assumption. This is important for 
understanding the extent to which the symmetry-based 
picture we develop here can be applied to realistic experi- 
ments. We then use several different numerical examples 
to illustrate how the resonant triad interactions most rel- 
evant to pattern formation may be controlled through a 
judicious choice of forcing parameters. 

The remainder of this paper is organized as follows. In 
Sectional we review basic ideas concerning the impor- 
tance of resonant triad interactions to Faraday wave pat- 
tern formation, including a discussion of some of the pre- 
vious theoretical and experimental work. Section lllll con- 
tains our symmetry-based analysis. We enumerate the 
most important weakly damped modes, calculate their ef- 
fect on pattern formation, and determine the dependence 
of this effect on the forcing parameters. Section Hvl con- 
tains a general discussion of the symmetry-based results. 
We study their range of validity with respect to 7 by 
comparing the symmetry-based predictions to numerical 
results obtained using the Zhang- Viiials Faraday wave 
equations 21]. In Section [3 we apply our symmetry- 
based results in several examples. In the first application, 
we consider weakly-nonlinear periodic patterns forced by 
a repeated sequence of (5-functions of alternating sign. 
In accordance with the results first reported in |22l |. we 
demonstrate how, by varying the spacing between the 
pulses, we may control the amplitude of the pattern. In 
the second application, we show how to construct a five- 
frequency forcing function which leads to dramatic sta- 
bilization of a complex pattern, namely an SL-I super- 
lattice pattern of the type observed in |7J . In the third 



example, we conjecture about a seven-frequency forcing 
function which should be conducive to the experimental 
observation (as yet, lacking) of 14- fold quasipatterns. We 
summarize and conclude in Section IVII 



II. BACKGROUND 

We lay the groundwork for our new results by repris- 
ing basic ideas from [H 113, 111 IH IH H^ on the role 
of resonant triads in Faraday wave pattern formation. 
We consider Faraday waves on an unbounded horizontal 
domain subjected to an arbitrary periodic forcing func- 
tion f{t). We use a dimensionless time t such that the 
common frequency is one, and expand /(r) in a Fourier 
series: 
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(1) 



where u = m,n,p, . . . are the forcing frequencies (dis- 
tinct and co-prime), |/„| are the forcing amplitudes, and 
0„ = arg(/„) are the corresponding phases. Without loss 
of generality, we take m to be the "dominant" frequency, 
i.e. we assume that fm (to lowest order) is the component 
that drives the critical modes (this does not necessarily 
mean that |/„i| is the largest of the \fu\)- There exists 
a bifurcation point |/,„| = |/m|'^'^** which depends on the 
physical properties of the fiuid, and on the other /„, be- 
low which the fiat fiuid state is stable to perturbations of 
all wave numbers, and at which perturbations of (generi- 
cally) one critical wave number kc become neutrally sta- 
ble. We consider the properties of resonant triads in a 
vicinity of this bifurcation in parameter space. 

Three wave, or triad, resonance is the simplest non- 
linear mechanism by which different waves may interact. 
The three waves involved have Fourier wave vectors kj, 
j = 1,2,3, satisfying 



ki 



(2) 



In this paper we are interested in the infiuence of the 
damped modes that are driven nonlinearly (through res- 
onant triad interaction) by the critical modes. Hence 
two of the wave vectors have the critical value |ki| — 
|k2| = kc. These waves, to first approximation, respond 
subharmonically to the dominant forcing component m 
and thus oscillate with predominant frequency m/2. The 
third wave in the triad has wave number jkaj = k^ and 
is associated with a damped mode with dominant fre- 
quency Q. The values of 51 most relevant to Faraday 
wave pattern formation arc determined in Section IIIII 
The condition ^ defines an angle of spatial resonance 
^ros G [0, tt) between the two critical modes: 
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We exclude the cases 0rcs = ""/S and ^j-os = 2tt/3 since 
these correspond to hexagons and k^ would then not be 



damped. Figure Q] shows Fourier space diagrams corre- 
sponding to the resonant triad described above. 



manifold reduction can be used to eliminate At. We find 
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FIG. 1: Fourier space diagram of spatially resonant triads sat- 
isfying ^ . The two neutrally stable modes have wave number 
|ki| — |k2| = fee and oscillate with dominant frequency m/2. 
The damped mode has Iks] = kd and oscillates with dominant 
frequency fl. (a) kd < kc- (b) kc < kd < 2kc. 

In the presence of damping, the primary Faraday insta- 
bility leads to standing waves (SW). We associate with 
the wave vectors k^ three complex amplitudes Aj which 
describe the slow-time evolution of the three standing 
wave modes pictured in Figure ^ (the fast-time subhar- 
monic oscillation of the waves has been factored out; see, 
e.g., iV^). Evolution equations for the Aj can be ob- 
tained by applying a standard reduction procedure to the 
governing equations, as we shall do in Section How- 
ever, at this point we are concerned only with the form 
of these equations, which is determined by the spatial 
symmetries. The action of spatial translation is 



Te : Aj ^ Aje''^ , 

e= (01,02), 01,02 e [0,2^), 03=01+02, 

while a reflection about ka leads to 
K : Ai <~i A2, 
and a rotation by tt acts as 



7^: A, 



A, 



(4) 



(5) 



(6) 



Equivariance under these three symmetries (see, £.17., 
[25J) requires that the differential equations describing 
the dynamics of the Aj take the form 



Ai = AiAi + aiA2A3 

+ {a\A,\' + bo\A2\^ + h\A3\')Ai, 

A2 = Ai^2 + ai^i^3 

+ ia\A2\^ + bo\Aif + bi\A3\^)Ai, 

A3 = A2A3 + a2AiA2 

+ (62|Ai|2+62|A2p+ 631^312)^3, 



(7a) 
(7b) 
(7c) 



to cubic order. The dot represents differentiation with 
respect to a slow time scale. All coefficients are real. 

Because Ai and A2 are neutrally stable modes and ^3 
is linearly damped (i.e., Ai = and A2 < 0), a center 



C^2 
A3 = -—A1A2 
A2 



(8) 



in a neighborhood of the origin. The (unfolded) bifurca- 
tion problem, to cubic order, becomes 

ii = AiAi+a|Ai|2Ai+&(0.,s)|A2pAi, (9a) 

A2 = AiA2+a|A2|M2+6(0re«)|AipA2, (9b) 



where 



6(0rcs) 



aia2 



(10) 



The coefficient 6(0) is the cross-coupling coefficient for 
SW oriented at an angle relative to each other and, 
above, it is evaluated at the angle of spatial resonance 
= 0rcs- The resonant contribution 6res arises from the 
presence of the damped kd mode. 

The resonant angle 0rcs ranges from to tt as kd varies 
from 2kc to 0. When kd is such that the natural fre- 
quency flikd) of the damped mode equals (or is nearly 
equal to) one of the special values that promotes a strong 
nonlinear coupling (as determined in Section lTTH) the con- 
tribution 6i-Gs to 5(0res) can be significant. This typically 
happens when ai and a2 become large in magnitude, 
and/or when A2 becomes small in magnitude. The res- 
onant contribution will then have a major effect on the 
stability of associated patterns. 

Consider further the system Q which has as steady- 
state solutions the trivial state \Ai\ — \A2\ ~ 0, the 
symmetry-related "striped" states |^i| > 0, \A2\ = and 
1^2! > 0,1^41 1 = 0, and the "rhombic" mixed-mode so- 
lution 1^1 1 = \A2\. We assume that a < 0, so that 
the bifurcation to the striped state is supercritical. A 
straightforward analysis yields the following stability re- 
sults summarized by Figure|2] For 6 sufhciently negative, 
i.e., 6 = bo + brcs < a, the (supercritical) branch of rhom- 
bic states with angle 0rcs is unstable. If b is increased 
such that |6| < \a\ (typically due to b^cs > balanc- 
ing 60 < 0) then the two modes mutually enhance each 
other's growth, and the rhombic pattern is stable. If b 
is increased further (due to an even larger, positive 6ics) 
such that b > —a, then the rhombic state bifurcates sub- 
critically. However, with the addition of fifth order terms 
(or higher) it is possible, even likely, that for the subcrit- 
ical case, the unstable mixed-mode branch turns around 
at a saddle-node bifurcation and creates a branch of sta- 
ble, finite amplitude rhombic states. Thus, we do not 
want to be unduly limited by the form of ||51). In the ini- 
tial stages of the pattern selection process, when modes 
on the critical circle are beginning to grow and compete, 
there will surely be an advantage for combinations that 
mutually enhance each other's growth. For these reasons 
we say that triad interactions contributing bros > are 
enhancing and those giving b^cs < are suppressing. 

The above example is just one very basic instance of 
the importance of resonant triads. In fact, triad reso- 
nances have implications far beyond the (in)stability of 




FIG. 2: Three qualitatively different phase portraits corre- 
sponding to ® with a < 0, Ai > 0. Top: b < a. Middle: 
a < & < —a. Bottom: b > — o. 



rhombic patterns. They may affect the stabihty of pat- 
terns within the framework of any Faraday wave bifur- 
cation problem possessing a subspace with d yna mics de- 
scribed by ®; see, for instance, pl.ll7lll9ll20ll2lll23|. 
In these cases, the logic is the same: &rcs > enhances 
patterns involving the resonant angle ^rcs while &rcs < 
suppresses them. 

The triad interactions discussed in this paper have im- 
plications for one-dimensional waves as well. In this case, 
with ki —li.2, the resonance condition Q becomes sim- 
ply 



ZKc — K(i , 



(11) 



which is the familiar 1 : 2 spatial interaction. When the 
natural frequencies of the two waves are such that a 
strong nonlinear coupling is allowed (as we detail in Sec- 
tion ^Oj we expect additional contributions to the cubic 
self-interaction coefhcient a in the SW equation 



dAi 

It 



AiAi +a\Ai\'^Ai 



(12) 



which is simply @ restricted to one spatial dimension. 
Since there is no spatial angle 9 to tune, we may arrange 
for a resonant situation (|ll|l by varying parameters in 
the dispersion relation, as in |19ll23l |. 



III. SYMMETRY CALCULATIONS 

We use the approach developed in [13, U^ l2J| to de- 
termine how the resonant contribution fej-cs to the cross 
coupling coefficient 6(0rcs) in I® depends on the forcing 
function Q . We consider a system of six traveling wave 
(TW) modes (see also [ly|) having the same wave vectors 
kj as the three SW modes described in Section [^ It is 
advantageous to consider TW first because the action of 
the temporal symmetries on TW is simple while on SW 
it is not. In this way, we make full use of the temporal 
symmetry and Hamiltonian structure before reducing the 
TW equations to the desired SW equations by means of 
a center manifold reduction. We thus expand the fluid 
surface height h(x,t), x G M^ as 



/i(x,i) = ^^Z±(t)e'(''^-"±^^-*) +C.C., (13) 

j=i ± 

where Z, are the slowly varying amplitudes and 

1371 = ^^2 = m/2, W3 = fl. (14) 

Spatial and temporal symmetries constrain the equations 
for the evolution of Z^ , as we now detail. 



A. Spatial symmetries 



Spatial translation symmetry acts on the TW ampli- 



tudes Zf as (c/. Eq.lU 



Te-.Z 



e= (01,02), 01,02 6 [0,2^), 03=01+02. 

A reflection about ka acts as (c/. Eq. [SI 



K : Zi ^-» Zn 



VI -^ ZV2 , 

and a rotation by tt induces ( c/. Eq. |^ 



7^: Z? 



ZJ- 



(15) 



(16) 



(17) 



We enforce equivariance under (|15|) - (|17|l to obtain the 
form of the TW amplitude equations to quadratic order. 
This truncation is sufficient to determine the leading or- 
der resonant contribution 6rcs to 6(0ros) in lO. We have 



Zi — L\Zy + L2Z^ + Q1Z2 Zg 



.3 (18a) 

+ Q2Z2 Z^ + Q3Z2 Z^ + Q4Z2 z^ , 
Z+ = L3Z+ + L4Z^ +Q5Z+Z+ (18b) 

+ Qe^i Z2 + QeZ^ Z2 + QjZ^ Z2 , 

where the remaining four equations follow from the ap- 
plication of (Cnil and ((T7|) . 

We now apply a standard reduction procedure to H18|) 
and compare this result with the SW equations l|njl. To 



facilitate the subsequent calculations we first introduce a 
phase shift to the amplitudes: 



■"1,2 






where 



ip = ip2-^i+-n, 



(19) 



(20) 



with ipi^2 denoting the phases of the coefficients Li and 



(i.e., Li 



|Li,2|e*'^i.= ). The TW equations ^ 



may be compactly written in the form 

Z = LZ + N(Z), (21) 

where Z = (2+ , Zf , Z+ , Z^ , Z3+ , Zg" )^ . 

The bifurcation to SW occurs when \L2 
we will see in the next section, IL2I ~ |/i 
bifurcation condition serves to define the critical am- 
plitude of the dominant forcing component \fm\- The 
critical eigenvectors are vi = (1,1,0,0,0,0)'^ and V2 — 
(0,0, 1, 1,0,0)"^. We use a multi-scale perturbation cal- 
culation to accomplish the reduction to SW, writing 



Li\. As 
ml, SO this 



I/. 



z = 


= T]{AiVi+ A2V2) +r]'^Z2 


m\ ~ 


- \frnr'+vV2\ + ---, 


d 
dt 


2 d 



(22a) 
(22b) 

(22c) 



where 77 ^ 1 is a small bookkeeping parameter and A12 
are the time-dependent SW amplitudes. At 0(77) the 
linear problem is recovered. At 0(7]^) Z2 is determined. 
At Oiji^) a solvability condition yields equations for the 
slow variation of the SW amplitudes: 



i2= AiA2+5res|AiPA2. 



(23a) 

(23b) 



The coefficients a and bo in © do not appear above be- 
cause the cubic order terms were omitted in ()18|l . For 
the purposes of this paper, we need only point out that 
the "nonresonant" coefficients a and bo are both 0(7) 
[12. I19I |2J| (recall that 7 is a dimensionless measure of 
the damping). The resonant contribution is given by 



f^rns 



Re{LiS} 
Rc{Li} ' 



(24) 



where 



S^QiWe~''^ + Q2We-''^ + Q3W + QiW, (25) 



with 



W - {LiU-L3U)/{\L3\^-\Li\^) 



U 



e"^Q5 + 2Qe 



'Qi. 



(26) 



and <p defined by (|20|l . Our analysis applies when jisl > 
IL4I, i.e. when the ka mode is linearly damped. 



B. Temporal symmetries 

Temporal symmetries constrain the coefficients 
Li, . . . , L4 and Q^, (. = 1, . . . , 7 in ifTH)) . In the absence 
of damping and forcing, the problem has an exact time 
translation symmetry 



Ta* : Z^ 



r^zk zLizuj At 



(27) 



with zuj given by l|14|l . and an exact time reversal sym- 
metry 



a : t ^ — i, 



Zl 



^7- 



(28) 



In the presence of finite damping and forcing, these tem- 
poral symmetries are broken. Nonetheless, they can be 
recast as unbroken parameter symmetries by allowing an 
appropriate transformation of the forcing parameters /„ 
and the damping 7. With this generalization the time 
translation symmetry 1)27(1 becomes 



TAt : Zf -^ z±e±'-^-^*, /„ ^ /„e 



■j,At 



and the time reversal symmetry (|28|) becomes 
a : {t, 7) ^ -(i, 7), Zf ^ Zj, U ^ /„ 



(29) 



(30) 



The damping and forcing are both assumed to be 
small, and are of the same order, i.e., |/„| ~ 7 <^ 1. 
A Taylor expansion of the coefficients Li, . . . , L4 and Q^ , 
consistent with H29() and (|30|l . leads to 



Li = -Vr-^, 


(31a) 


L2 = -iKfm, 


(31b) 


L3 = -Qr7, 


(31c) 


Li = -~iy.iF2n, 


(31d) 


Qi = iqiFi, 


(31e) 



where only the leading order terms have been kept. The 
expansion coefficients are all real, and Vr, Qr > since 
they correspond to damping terms. The factor of /„ 
in the expansion of L2 reflects the fact that the critical 
modes respond parametrically to the dominant compo- 
nent /„. The factor F2Q. in the expansion of L^ rep- 
resents an analogous parametric forcing term for the 
damped mode composed of products of the /„ and fu 
whose frequencies sum to2Vt. When 20, forcing is present 
in 1^, then, to lowest order, F20 = /2n; otherwise L4 = 
at 0(7). 

The Fi in (|31e|) . in accordance with (|29|l . must con- 
tain products of the /„ (and /„) whose frequencies 
are such that (Qi, Q5)e'^*^"™'*, (Q2, Q7)e-*(™+"^*, and 
[QciiQiiQo)^ are time translation-invariant quanti- 
ties. Since we are interested in understanding when the 
effect of resonant triads is significant, we focus on those 
cases where &rcs is 0(7) or larger; this requires that one 
or more of the quadratic coefficients Qi is 0(7) or larger. 
A straightforward calculation shows that this scaling can 
occur only if fJ G {m, 2m, n,m ± m,n — m}, Q > 0, for 



some frequency n, and we henceforth restrict attention 
to these cases. Note that, since there are many frequen- 
cies in /(r), these sets of relevant D, values can overlap. 
For instance, an "m — n mode" is also a ^^p — m mode" 
ii n + p = 2m, n < m. An important (and somewhat 
surprising) result of our symmetry calculation is that the 
contribution 5rcs arising from a given damped mode with 
frequency Q involves (at leading order) a maximum of five 
frequencies: the dominant frequency m, up to three other 
frequencies appearing at 0(7) in the three sets of cou- 
plings coefficients {Qi,Q5), iQ2,Q7), and {Q3,Qi,Q6), 
and potentially one more frequency, 2f2, that parametri- 
cally forces the damped mode at 0(7), thus making L4 
nonzero at leading order. The effect on the triad interac- 
tion of any additional forcing components in |^ will be 
0(7^) or smaller. 



C. Results 



We combine the results of Sections IIII Al and IIIIBI to 
obtain leading order expressions for foj-cs in with ex- 
plicit dependence on the damping 7, the forcing ampli- 
tudes \fu\ and the forcing phases 0„. For each fl there 
are a finite number of qualitatively distinct cases to con- 
sider. These are distinguished by the number of relevant 
frequencies involved (up to five) and the manner in which 
they enter the problem (through L4 and the Qe). Hav- 
ing chosen one of these, we substitute the corresponding 
expressions for Qi, and the expansions for the remain- 
ing T W coefficients shown in (I31II , into equation (|24|l for 
6rcs- The results are summarized in Tabled and will be 
discussed in the next section. 

To make this table of results manageable, we make use 
of the following definitions: 



as 



91 95, 

9297, 

295(93 -94), 
9197-9295, 



(32a) 
(32b) 
(32c) 
(32d) 



{29196 + 95(93 -94)} Ai/|Aj|, (32e) 

{29296 -97(93- 94)} A,/|A,|, (32f) 



and 



P20($) = 



1^3! +M»l/2n|sin$ 
|i3P-|M./2nP ' 



(33a) 



R (<b <b \- l-^3|sm$i-H^,|/2n|cos^2 , , . 

^2o(<±'l,<±'2) - T-j—r^ \ 7 — [^ • {■^■^°) 

In the above, the qj and Xi are defined by (|31|) . The 
relevant phases <&, $1, $2 appear in Table Q] 



D. Hamiltonian structure 

We now discuss the implications of Hamiltonian struc- 
ture in the undamped problem (see IH HI mUl IH 



ISOl W\). This is a stronger assumption than that of 
time reversal symmetry l)29|l alone. We suppose, as in 
El 113, 111, that the undamped TW equations ^ can 
be derived from a Hamiltonian Ti.. Because the ampli- 
tudes Z and Z need not themselves be canonically 
conjugate Hamiltonian variables, we write Hamilton's 
equations in the generalized form 



^t-T 



1 dn 



rjdZ, 



T' 



n = r2, r. G 



(34) 



This takes account of scaling transformations like 
{Z^,Z^) -^ ri{Z^,Z^), Zf -^ raZ^ that preserve the 
Hamiltonian character of the dynamics, and are needed 
to relate the underlying canonical variables to Z and Zj^ 
in H18|l . For inviscid Faraday waves the surface height 
h and the surface velocity potential are the underlying 
canonical variables (see, e.g., |23, 123). Using this fact 
we find that, to leading order, rf = r2 = m/{2kc) and 
r| = Vt/kd are appropriate prefactors in P^ (see |23| 
where a similar factor arises in the corresponding canon- 
ical transformation). 

Requiring that 7i be a real-valued function, invariant 
under the symmetries ((T3|) - ((T7|l . if^ and (|Sn|l . we find 
that the equations of motion (|34|l are equivalent to H18|l 
only if qi = rq^, 92 = ?'97, and gg = 94 = rge with 
r = r'^/rf. These conditions imply, for the results in 
Table n that 



ai > 0, a2 > 0, as = 0, a^ = 0. 



IV. DISCUSSION 



(35) 



We now discuss TableQlin some detail, highlighting the 
most important features of the results collected there. 
We then investigate the range of validity of these results, 
which were derived under the assumption of weak damp- 
ing 7. To do this, we introduce the Zhang- Vifials Faraday 
wave equations and use them to perform explicit numeri- 
cal calculations that demonstrate the range of 7 for which 
the symmetry-based results provide an accurate predic- 
tion. 



A. Highlights of results 

Some general comments on the organization of Tabled 
are in order. Note first that there are many cases which 
do not need to be listed because they can be obtained 
simply by relabeling the different frequencies. For exam- 
ple, the case {m, n,p, q,- ;fl) = (m, n, 2m+n, m+n, ■ ; ra-\- 
n) is equivalent to the case (fourth up from the bottom 
in Table P) {m, n, p,q, ■ ;il) — (to, n,m + n,n — m, ■ ;n) 
with n -f-^ q. 

There are six groupings in the table. The first shows 
the five important damped modes and their contribution 
to b-cca when there is only one type of coupling at 0(7) or 



TABLE I: Leading resonant contribution 6ros to b in Q for the most important damped modes. For a damped mode with 
frequency f2, there are at most five forcing frequencies m,n,p,q,r which affect &rcs- Here, Tn,n,p,q,r,fl > and x G Z^. Each 
expression for (m, n,p, q, r), given Q, is excluded from those of entries further down the table, in which additional relationships 
hold. Dots indicate an arbitrary commensurate frequency, if present, which does not affect fores at lowest order. Entries whose 
listed frequencies have a common factor {e.g., those with x) are assumed to be part of a forcing function with other, relatively 
prime, frequencies. For * the ± follows sign(m — n). See (I32II and 13311 for definitions of qi, . . . ,a6, P2n, and R2n used below. 
Certain entries are reproduced from I2Q|; the cases that involve more than three forcing frequencies are new. 



{m,n,p,q,r) 



n 



(m, ■, 
(m, ■, 
(m, n, 
(m, n, 
{m,n, 



■) 



m 

2m 

n 

m iin 

n ~ m 



{m, 2m, ■, ■, ■) m 

(m,Am, ■, ■, ■) 2m 

(m,n, 2n, ■, •) n 

{m,,n,2m ±2n,-,-) m,±n 

(m,n,2n — 2m,-,-) n — m 



{m, 2m, -, ■, ■) 2m 

{m,3m, -, ■, ■) 2m 

[m,n, \m — n\, -, -) n 

{m,n,m + n, -, -) n 

{m,n,2m±n,-,-) m±n 



(3a;, 2a;, 4a;, ■, ■) x 

[m, 2m,, Am,, -, •) 2m 

(m, 3m,, 4m, •, •) 2m 

{m,,n,\m ~ n\,2n,-) n 

(m, n, 771 + n, 2n, ■) n 

(m, n, 2m ± n, 2m ± 2n, •) m±n 



(m, 2m, 3m, •, ■) 

(m, n,m + n,\m — n\, ■) 



2m 



(3a;, a;, 4a;, 2x, •) 



(m, 2m, 3m, Am, •) 



Leading resonant contribution bi. 



— ai/|L3 
-ai|/,„P/|i3 
— aal/nl 71-^3 

-ai|,fn|Vl^3| 

a2|/n| /\La\ 



-QlPn($) 
-Qll/„|'Pn($) 

-a3|/„pPp($) 
-ai|/npPp($) 

«2|/„pPp($) 



(-ai|/m|^ -a3|/n|''+Q5|/m||/n|sin'l>)/|L3 

(-ai|/mP + a2|M'+a4|/™||/n|cos$)/|L3| 

(-ai|/pp-a3|/„|2+a5|/„||/^|sin<l.)/|L3| 

/•„ I -f |2 „ 1 -f |2 I „ I i- II i- I „;„ /RN /I r 



(-ai|/pp-a3|/„|2+a5|/„||/^|sin<l.)/|7 
{a2\fpf ~as\f„\' +ae\fn\\U\sm$)/\L, 
{ailfpf -ai\fn\^ + a4\fn\\fp\cos$)/\L3 



-ai|/p|^Pp($i-$2)-a3i/„l"Pp($i+<E>2) 

+ «5|/n||/p|i?p($l,$2) 
-ai|/npPn(3'l+$2)+a2|/ppPri($2-$l) 

+a4|/„||/p|Pn($i- 90°, $2+90°) 

-ail/,„pPp($l- $2) - Q3|/nPPp($l+ $2) 
+ Q5|/^||/n|Pp($l,$2) 

-ai|/,„|'Pp($i- $2) + a2|/„pPp($i+ -1>2+ 180° 

+ a4|/m|l/n|ii:p($i+ 90°, $2+ 90°) 

-ai|/p|2p,(2$i- $2) - a3|/npP,($2) 

+ «5|,/n||/p|P,($l,<E>2-$l) 
a2|/ppP,(2<E>l-$2)-Q3|M'P,('&2) 

+ a6|/n||/p|P,($l,$l-$2) 

-ai|/„|2p,($2-2$l)+Q2|/p|'Pg($2) 

+ a4|/n||/p|P,($l+90°,$2-$l-90° 



{(a2|,fp|^ - Ql|,fm|^ - Cl3|/n|^ + Q4 |, fm | | ,fp | COS $1 

-I- (Tc I /■„ 1 1 /■„ I sin <i)n 4- (Tel f_ II /■_! sin ('(60— <i)> H/l /",o 



|/n||/p|sin($2-''l>l)}/|L3| 
2 . . 



XV'-^2|./p| «-ti|./m| — *-^ 

+ a5|/ml|/n|sin'l>2 + Q6|jn||jp|-"n 

{(a2|/pP - ai|/,P ~ a3|/nP + a4|,/-p||./-,, .™ .. 
+ ct 5\.U\\fq\sm{^i+ ^2) + Q6 1. fn 1 1, fp I sin $2}/ 1 1/3 

rTTTTTT^ a; N i " 



I COS<&i 



relevant phase(s) 



$ = 0n— 2(/!>m 

$ = 3<^„ - 2,^„ 
$ = 2<^„ — 0p 



$ = 



2(j)mT 2(/>n 



<& = <^p+ 20m— 20„ 



$ = 071— 20„ 
$ — 4>n— 30„ 



$ = 
$ = 

$ = 20„ 



- 0m ± I 
I + 0n — I 



$1 
$2 
$1 
$2 

$1 
$2 

$1 
$2 

$1 
$2 

*1 
$2 
$1 
$2 



4>n— 4'm+ 4>p 

4>m+ <t>n~ 20P 

4>n + 0P — 20m 

20n — 0p 

0,1 — 20m 

0„ + 20m — 0p 

0n — 30m 

0m + 0n — 0p 

0n — 0m ± 0p * 

20„— 0g 

0m + 071 — 0p 

2071— 09 

20m ± 071 — 0p 

20m— 20P+ 0g 



2m 



(m,n, m + n, |m — n|, 2n) n 



-Ql|/,l"P,($2-$l)+a2|/p|"P,($l+'l>2) 

-a3|/7i|'P,(2$3"$i-'l>2)+a4|/p||/,|P,($i+90°,$2-90°) 

+ a5|/7i||/5|P5($3-$l,$3-$2)+a6|/7i||/p|P,($3,$l+$2-$3) 
-ai|/mpP,(--l'l-$2)+a2|/p|'P5($2-$l) 

-a3|/7i['P,(2cE.3-$i-'l>2)+a4|/m||/p|P,($2+90°,$i+90°) 

+ «5|/7n[|/7i|P,($3-$l-$2,$3)+a6|/7i||/p|Pg(<E>3-$l,4>3--l>2; 
-Ql|/,pP.($2-$l)+a2|/ppP7-(<E>l+$2) 

-«3|M'P7-(2$3-$i-$2)+a4|/p||/,|P7.($i+90°,$2-90°) 

+ "5|/7i||/q|Pr($3- ■I>l,$2-$3) +Q6|/7i||/p|Pr('I>3,'I'2+'I>3- $l) 



$1 
$2 

$1 
$2 



Ap — 30m 

Ati — 20m 

i>p — 20m ± 09 * 

Am + 071 — 0p 



*1 
$2 
$3 

$1 
$2 
$3 

$1 
$2 
$3 



: 20m — 0p — 0q 
20g — 0p 
0m + 071 — 0p 
0m + 0p — 0IJ 
30m + 0p 
20m + 071— 0q 
20m— <j)pT 4>q* 
0r— 0p± 0IJ * 
0m + 071 — 0p 



lower and no parametric forcing /2f2- In these cases there 
is no (leading order) dependence on the forcing phases 
^„. In the second section the same damped modes have 
been parametrically forced. The factor l/lL^l is then 
replaced by P2n(*&) of |j32il)- This is a strictly positive 
oscillatory function (I-L3I > |/ii/2n| for damped modes) 



with extrema at $ = ±90°. The third and fourth sec- 
tions are analogous to the first and second, but with two 
types of coupling rather than one - similarly for the fifth 
and sixth sections, but with all three possible quadratic 
couplings {i.e., all Qe are linear in the /„). 

Two of the damped modes appearing in the table war- 



rant special mention. The 17 = m mode stands out be- 
cause its influence is especially strong. For this mode, the 
largest quadratic terms in l(TH|l are 0(1), and the result- 
ing contribution 6rcs is 0(7""'^). In contrast, for all of the 
other damped modes, the strongest quadratic couplings 
take place at 0(7) and lead to 5res of 0(7); these 0(7) 
contributions are of the same order as a and &o in lO, 
but can still have significant effects on pattern seelction, 
as demonstrated in Section FVl 

The second special case is the Vt = 2m mode. Although 
this mode satisfies all the necessary temporal constraints 
to make a significant contribution fores, it cannot enter 
into resonant triad interactions with the critical modes 
because its wave number is too large, i.e., k^i > 2fcj, and 
(PI cannot be satisfied; one can estimate the relevant 
wave numbers from the inviscid fiuid dispersion relation 
(see 



|19j and Section [IV B|l . However, this mode may 
have relevance for other systems such as ferrofiuids in a 
magnetic field where the dispersion relation is nonmono- 
tonic |33j, and hence we have kept it in the table. 

A key result of Table Q] is the important role played by 
the relative phases (/)„ in the forcing function Q . For all 
but the most simple cases (in the first section of the ta- 
ble) , b-cca depends on combinations of the forcing phases 
which are invariant under the time translation symmetry 
Tr of (P5|l : one phase is always arbitrary, associated with 
the choice of origin in time, while any physically mean- 
ingful phase must be invariant under (|29|l . This phase 
dependence provides a very convenient way to tune the 
strength of the nonlinear interactions, as the numerical 
examples of Section M will demonstrate. 

Another important aspect of Table ^ pertains to the 
sign of 6res- Recall from the discussion of Section ITTI that 
if 6res > interactions involving critical modes separated 
by the angle ^i-cs will be enhanced, whereas if 6ros < 
they will be suppressed. Relations (|35|l mean that for 
simple couplings (the first two sections of Table P) the 
sign of 61CS is determined, and thus one knows which ef- 
fect (if any) to expect. In particular, the ^ — m,^ = 2m, 
and U, = m ±n modes are suppressing while the D, = n 
mode is inconsequential. The U, = n — m mode, in con- 
trast, is enhancing, and thus is of great interest because 
it may be used directly as a selection mechanism. The 
effect of this difference frequency mode on pattern selec- 
tion was examined in [13, and indeed, it is likely respon- 
sible for stabilizing the superlattice pattern observed in 
\j\. We examine the difference frequency mode further 
in Section Ivl 

A final noteworthy feature of Table|l]concerns the effect 
of parametrically forcing the damped mode with a fre- 
quency 251. A comparison of the factors l/jLaj, P2o($), 
and i?2o(*I'i, $2) reveals the potential for a small denom- 
inator in the latter two cases. The parametric forcing can 
increase |6res| and amplify the effect of the damped mode 
provided this denominator does not become excessively 
small, which would indicate that the damped mode is 
nearly critical and that the reduction leading to H23|) is 
breaking down. This feature will be exploited as well in 



some of the examples of Section El 



B. Zhang- Vinals hydrodynamic equations 

In this subsection, we investigate the range of damp- 
ing 7 for which our symmetry-based results are valid. To 
carry out this investigation we perform explicit numer- 
ical calculations using the Zhang- Vinals hydrodynamic 
equations (introduced below). In particular, we use the 
method described in [13 to calculate the cross-coupling 
coefficient 5 in © as a function of 0, the angle between 
ki and k2 in Figure Q] It is sufficient to take 9 £ [0,90) 
since 5(180° -6)^ 6(6*). 

The Zhang- Vinals equations |22 describe the dynam- 
ics of small amplitude Faraday waves on a deep, nearly 
inviscid fiuid layer. We use the same scaling of the equa- 
tions as in 23j, writing them in the form 

{dr ~ "fV^)h - V <P = J'(/i,$),(36a) 
{dr - 7V2)$ - (FoV^ - Gir)) h = Gih, $), (36b) 

where G{t) = Gq — /(t) and the nonlinear terms are 
given by 



+ V i hV{hV<i>) + -h^V^<^ 

gih,<^) = l(5<i>)2-i(v$)^ 



(37a) 



(37b) 



-^ToV-{iVh)iVhf}. 

Here /i(x,i) is the fiuid surface height, ^(x, t) is the sur- 
face velocity potential, and x is the two-dimensional spa- 
tial coordinate. The operator V multiplies each Fourier 
component of a field by the modulus of its wave number. 
I.e., Ve'^'' = |k|e*'^''. 

The equations depend on three dimensionless fiuid pa- 
rameters: the damping parameter 7, the gravity number 
Go, and the capillarity number Fq. These fiuid parame- 
ters, and the dimensionless forcing amplitudes /„ in (Q 
are related to the physical parameters by 



_ 2iyk'^ _ gok 

1 = , Go = -^, 



ak^ g^k 

To = — ni In = ^-5-- (38) 



Here v is the kinematic viscosity, a is the surface ten- 
sion, p is the density, and to and the gu are the Fourier 
amplitudes in the original (dimensioned) forcing function 



5W = E 



iuujt I 

gii,e + c.c. 



gu e 



(39) 



Additionally, k satisfies the inviscid gravity-capillary 
wave dispersion relation 



~ ak^ fmuj 



(40) 



and (7o is the usual gravitational acceleration. Note that 
Go and To are not independent parameters since H.S8(1 and 
(gni imply that 



Go + To — 



(41) 



The dimensionless dispersion relation (c/. Eq. I40|l also 
gives the natural frequency ^(k) of undamped, unforced 
waves as a function of their wave number k: 



n^ = Gak + Tok^ 



(42) 



For small damping 7, (|42|l provides an excellent esti- 
mate of the wave number which oscillates at a given fre- 
quency, even for forced waves; we make use of this fact in 
Section Since the critical modes oscillate with dom- 
inant frequency to/2, we have kc ~ k{m/2) = 1, where 
fc(ri) is the inverse of the dispersion relation from (|42|l . 
One may then choose a damped mode with frequency Q, 
find fc(il), and then apply Q to estimate ^ics- 



C. Validity of symmetry-based results 

To investigate the applicability of our results for finite 
values of 7, we focus on an example using three- frequency 
{m, n,p) = (8, 7, 2) forcing and quantify the effect of the 
51 = 8 — 7 = 1 damped mode; this corresponds to the 
penultimate entry in the second section of Table Q] Al- 
though this mode does not necessarily lead to the most 
significant resonance, we study it as an instructive ex- 
ample to address general questions about the validity of 
our symmetry results. Damped modes which play a more 
important role are examined in the applications in Sec- 
tion |V| 

From the Hamiltonian considerations in Section fllll we 
have ai > 0, and thus b^-cB < 0. We set Fo — 16 in (|SS|l . 
fix the ratios of the forcing amplitudes at |/n|/|/m| — 0.4, 
\fp\/\fm\ — 0.08, and compute the coupling coefhcient 
b{9) using the method described in [l3. As predicted on 
the grounds of symmetry arguments, there is a dip in the 
plot of b{9) around the angle 9rcs ~ 23° where the fi = 1 
mode is in spatial resonance. An example is shown in 
Figure|21for 7 — 0.1 with $ (which appears in the fourth 
column of Table QJ set to 0. 

In the discussion that follows, we study various prop- 
erties of 6rcs as the damping parameter 7 is varied. In 
this discussion, it is important to realize that the results 
will depend on the chosen value of m, on which 7 in 
()38|l depends indirectly through (|40|l . When generaliz- 
ing the results shown below to other forcing functions it 
is, in fact, better to look at the quantity 7/771 (c/. I38|l . 



-3.04 



-3.07 



-3.10 



-3.13 




FIG. 3: Coupling coefficient b{9) in @ computed from 13611 
using three-frequency {m,n,p) = (8,7,2) forcing, 7 = 0.1, 

Fo = 16, \fu\/\fm\ = 0.4, \U\/\fm\ = 0.08, <Pm ^ (l>n ^ <Pp = 

0. 



This alternative nondimensional measure of the damp- 
ing utilizes the critical wavenumber and the dominant 
frequency (tow, as opposed to to) and is therefore better 
suited for quantitative comparison across forcing func- 
tions with very different to. values. We have used the 
scaling (|38l) . which utilizes the common frequency, to be 
consistent with previous work 19', 2(l l23Ll24| . 

We first consider the scaling of | ^ics | as 7 is varied with 
$ = 0. It follows from the result in Table |l] that 



f^rcs C5C I Jji 



\p4Pp 



(43) 



Furthermore, recall from H31I) that IL1.3I oc 7 and IL2I fx 
/,„. Since, at the onset of SW, \Li\ — \L2\ (see Sec- 
tion HTl, we have |/m| oc 7. Since |/m|, |/„|, and |/p| 
are held in a constant ratio, we also have |/„|, |/p| ex 7. 
Thus (|43|l becomes simply 6,cs oc 7. This scaling is con- 
firmed by the numerical results of Figure0]where we hold 
<& = and plot |5ies| as a function of 7. The numerical 
data is shown as points. For comparison, a line of slope 
one is drawn through the first data point, confirming the 
proportionality to 7. The theoretically predicted scaling 
holds reasonably well up to 7 ^^ 0(10^^), and the numer- 
ical result does not strongly diverge from the prediction 
until 7 « 0.5. 

Next, we examine the scaling of the half-width ^ of the 
dip dX 9 ^ ^rcs- For 9 ^ 9 res, the natural frequency of kd 
will differ from the resonant frequency (m, 2m^ n, etc.). 
At leading order, this detuning appears in the coefficient 
L3 as an imaginary part, i.e., L3 = —Qr^ + iQi {cf. Eq. 
131b ). If the detuning is small, the linear approximations 
Qi (X kd — krcs (X 6 — Ores cau be used (here fcres is the wave 
number associated with the resonant frequency) and so 
Qi — c{9 — ^res) for somc real constant c, i.e., 



L, 



-Qr"f + ic{9 



s). 



(44) 



10 





FIG. 4: Resonant contribution brca as a function of the damp- 
ing parameter 7. The dots correspond to a numerical com- 
putation using 1361 1. The straight hue of slope one confirms 
the fores oc 7 scaling predicted by symmetry arguments. The 
capillarity and forcing parameters used are the same as those 
in Figure El 



FIG. 6: The values of $ at which fo(6'rcs) takes on its minimum 
and maximum values as a function of the damping 7. The dots 
correspond to numerical data, while the lines at 90° and —90° 
show the predicted minimum and maximum respectively. The 
capillarity and forcing amplitudes used are the same as those 
in Figure 121 



Substituting this expression into the result from Table HI 
shows that ^ c>c 7. Numerical results are displayed as 
points on the log-log plot in Figured For comparison, we 
plot a line of slope one fit through the first data point. As 




FIG. 5: Half- width '^ of the resonant "dip" as a function of 
the damping 7. The dots correspond to a numerical compu- 
tation using l|36|l . The straight line of slope one confirms the 
predicted 4" oc 7 scaling. The capillarity and forcing parame- 
ters used are the same as those in Figure O 

with the dip magnitude |6rcs|, the theoretical prediction 
remains reasonable up to 7 ~ 0(10"^). 

Finally, we consider the dependence of b{6rcs) on $ = 
02 — 208 -|- 207, and examine how this <& dependence 
changes with increasing 7. From Table HI we expect 
that the dependence is sinusoidal and, from the fact 
that fii > |2J| for (|35|l . we anticipate 6(0rcs) reach- 



ing a maximum (i.e., having the shallowest dip) near 
$ = $niax = —90° and reaching a minimum (i.e., having 
the deepest dip) near $ = <l>niin = 90°. Figure shows 
how the numerically calculated values (dots) of $niin and 
*&max differ from the theoretical predictions (lines) as 7 
is increased. To elucidate the departure from the theo- 
retical prediction, we show three profiles corresponding 
to three different values of 7 in Figure [7| In Figure [71i, 
7 = 0.04 and the profile, as predicted, appears sinusoidal. 
The vertical lines drawn at ±90° agree reasonably well 
with the numerically calculated maximum and minimum. 
In Figure[7lD, 7 = 0.2 and, although the profile is still si- 
nusoidal, it is shifted by approximately 45° with respect 
to the theoretical prediction. In Figure [3;, 7 = 1, and 
the profile no longer resembles a sine function. This is 
demonstrated further by the plot in Figure \7\i, which 
shows the Fourier transform of the data in Figure [7|;. 
The zero component (i.e., the ^-independent part) has 
been removed, and the remaining data has been normal- 
ized so that the strongest component has magnitude one. 
The data indicate that higher harmonics of $ are now 
important. Note that the phase shift appears well before 
the higher harmonics come into play (see Figure Et"); a 
fact that can be understood as follows. The <& depen- 
dence in Tabled originates with the phase of terms in the 
normal form reduction, and depends on products of the 
coefficients in (|31ll . If the next order terms in the expan- 
sions describing these coefficients are kept, a phase shift 
of 0(7) is obtained. In contrast, higher harmonics of $ 
are generally associated with higher order (as opposed to 
next order) terms in the expansions (|31|l . This is a result 
of time-translation symmetry, which requires that terms 
involving additional powers of the forcing amplitudes /„ 
only appear in certain combinations. The specific order 
in 7 at which these new terms become relevant depends 
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FIG. 7: (a - c) Dependence of b(6rcs) on the phase $. (a) 
With damping 7 — 0.04. As predicted by the symmetry ar- 
guments in Section IIIII the phase dependence is sinusoidal 
with minimum and maximum near ±90°. (b) 7 = 0.2. The 
phase dependence is sinusoidal, but there is a phase shift of 
approximately 45°. (c) 7 = 0.1. The dependence is no longer 
sinusoidal, (d) Fourier transform of the data in (c). The 
zero component has been removed and the remaining data 
has been normalized so that the strongest component has 
magnitude one. The dependence on higher harmonics, e.g. 
2$, 3"1>, 4<1> is apparent. For all plots, the capillarity and forc- 
ing amplitudes are as those in Figure 



in complicated fashion on the particular choice of forcing 
frequencies. 

In this section we have explored the validity of our 
symmetry results with respect to the small 7 assumption 
under which they were derived. For small 7, the symme- 
try results are in excellent agreement with the numerical 
ones. For larger 7 , the scalings predicted by symmetry 
are not correct. However, many of the important qualita- 
tive features are preserved. In particular, even at larger 
7, increasing 7 increases |5res|- Furthermore, even though 
the dependence of &rcs on $ is no longer sinusoidal, there 
are still special phases $min and $niax which minimize 
and maximize feres; suggesting that even in experiments 
with large damping, tuning the forcing phases may be an 
effective means by which to control resonant triad inter- 
actions important to pattern formation. 



V. APPLICATIONS 

The results in Table|l]may be used to understand - and 
control - certain phenomena in Faraday systems. For 
each of the following examples, we apply our symmetry- 
based methods and demonstrate the results via numerical 
calculations using (|36|l . 



A. 1:2 temporal resonance and impulsively forced 
Faraday waves 

We focus on the cases for which H. = m in Table U 
so that the the critical modes and the damped mode 
are in a 1 : 2 temporal resonance. From the Hamiltonian 
considerations in Section IIIII ai > and thus feics < 
0. Also, recall from Section IIVI that for this case, the 
modes are coupled at 0{1). Therefore, the contribution 
6rcs is 0(7""'^), which is larger than for the other cases, 
where 5res is only 0(7). In short, the f2 = to mode 
has a very strong influence on b{6). The implications of 
this well-known resonance for Faraday waves have been 
investigated in a number of studies, including ,21.] . 

When /2m forcing is present, the size of 6res depends on 
the phase $ — (j)2m — '2-4>m] see the first entry in the sec- 
ond section of Table Q] This phase dependence has previ- 
ously been calculated in (l4| by means of a perturbation 
expansion on the Zhang- Vinals model H3t)|l . Our work 
confirms the phase dependence in a model-independent 
manner, exclusively by means of symmetry considera- 
tions. The phase dependence gives us a convenient and 
powerful means by which to control the 1:2 resonance and 
influence the shape of b{9). In particular, using $ = 90° 
maximizes the effect of the resonance, while $ = —90° 
minimizes it. 

In Figure IHl we show a numerical example for (to,, n) = 
(1, 2) forcing. The parameters in (|36|l are 7 — 0.008 and 
Fo = 0.125. The forcing amplitude ratio is |/n|/|/m| = 
0.396, which is far from the codimension-two point 
|/n|/|/m| — 3.53 at which waves with dominant fre- 
quency n/2 set in. The Q. = m mode has wave number 
k{m) « 1.83, and thus O^cs ~ 47.1°. Consistent with Ta- 
ble HJ a dip in h{ff) is found at this angle. As predicted, 
by choosing $ = 90°, we achieve the largest dip at 6'ios 
and thus a strong suppression of patterns involving an- 
gles near this one. On the other hand, using $ near 
—90°, actually reduces the effect of the triad interaction 
by a factor 1/2 < |L3|/(|L3| -I- \jlif2n\) < 1 relative to the 
single-frequency case, so the suppression is much weaker. 

As discussed in Section ^ the spatiotemporal reso- 
nances we consider in this paper may also affect the self- 
interaction coefficient a in the one-dimensional analogue 
of ©, namely (|12|l . In the case of the 1 : 2 temporal res- 
onance, the condition fl{kd) ~ 2ilkc must be satisfied 
along with (|ll|l . There will then be a contribution to 
the self-interaction coefficient a in H12I) whose dependence 
on the forcing and damping parameters is precisely that 
given in TablcQ] In practice, one may vary the frequency 
f2 by tuning the capillarity number Fq which appears in 
the dispersion relation (|42|l . In an experiment, this might 
be achieved by varying the base forcing frequency uj (see 
Eq.EHJ- 

The results of Table ^ for the 1 : 2 spatiotemporal 
resonance and its effects on the self-interaction coeffi- 
cient a may be used to understand certain features of 
impulsively- forced Faraday waves, i.e., waves forced by 
a periodic sequence of impulses rather than a smooth 
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FIG. 8: Effect of relative forcing phase on the first harmonic resonance, i.e., resonance with the Q — m mode, for {m, n) = (1, 2) 
forcing. The relevant phase $ is given in Table U (a) Cross-coupling coefficient b{0) with "1> = 90° and $ = —90°; the single 
frequency case (dashed line) is shown for reference, (b) Dip magnitude &(0ros) versus $. For these calculations, the parameters 
in Il3fcit are 7 — 0.008 and Fq = 0.125, and the forcing amplitude ratio is j/n|/l,/m| ~ 0.396. 
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FIG. 9: Schematic representation of the asymmetric 5- 
f unction forcing specified by H45|l . 



ares becomes more negative) as the asymmetry parameter 
a is varied across the interval (0, 27r). This observation is 
consistent with the results in Table Q] as we now explain. 
From Tabled there are at most two forcing frequencies 
which affect the D, = m damped mode at leading order, 
namely m and 2m. We therefore consider a drastic trun- 
cation of the Fourier series for the forcing function H45|l . 
keeping the first two terms, which are the only terms 
affecting the resonance at leading order: 



fit) = fie'' + /2P/^* + c.a 



(46) 



where 



forcing function of the form of Q . Impulsive forcing was 
studied first in 34] and subsequently in J23 • The original 
motivation for studying impulsive forcing is that it allows 
for a purely analytic linear stability calculation for arbi- 
trarily large damping, in contrast to the case of smooth 
periodic forcing functions for which the linear analysis 
must be performed numerically |l3l |35j and analytic ap- 
proximations are only valid for small damping. 
In |22|| the forcing function takes the form 



fit) = fsJ2 ^(* " 27r?i) - Sit - 2Tm - a), 



(45) 



n=0 



representing an alternating sequence of ^-functions of 
strength fs. The sequence has a temporal asymmetry 
controlled by the parameter a G (0, 27r), which deter- 
mines the amount of time between a positive pulse and 
the subsequent negative pulse. A depiction of l|15|) is 
shown in Figure In [22|, a(ro) is calculated from (|SS|l . 
and a large dip at Fq = Fj-es is observed, where Fj-cs is the 
parameter value for which the 1 : 2 spatiotemporal reso- 
nance is satisfied. For 7 sufficiently small, it is noted that 
this dip becomes more negative (i.e., the corresponding 
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(47) 



From (|46|l . we see that varying a affects both the am- 
plitudes and the phases of the forcing components. For 
(|36|l with weak damping and forcing, and for the two- 
frequency truncation (|46|) , the Faraday instability occurs 
when 1/1 1 = 7 (this follows directly from the results in 
[l9|). By setting fs equal to its critical value and making 
a translation in time, we can write the forcing function 
at onset as 



where 



/="*(i)=7e^*-f i^2e''*-H 



F2 = -2i7Cosf-j . 



(48) 



(49) 



The first entry in the second section of Table ^ indicates 
that the 1 : 2 spatiotemporal resonance produces a neg- 
ative contribution Oics to the self-interaction coefficient 
given by — aiP2($) where $ = arg(F2). Using the ex- 
pression H49|) and simplifying reveals that 



Ore 



Qfl 



\L3\ + 2\fi,\-fcosia/2)' 



(50) 
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which decreases as a is increased across the interval 
(0,27r). This is consistent with the observation in |2^ . 
which successfully fits numerical results to this functional 
form, at least for small 7. From (|12|l we see that the 
periodic striped state has a steady state amplitude of 
1^1 P = —A/a. Thus, experimentally, the wave height 
may be controlled by varying a. Larger a causes smaller 
a and, consequently, larger amplitude waves. 



B. Stabilization of superlattice patterns with 
multi-frequency forcing 

We now generalize the simple one-dimensional exam- 
ple just presented. Our symmetry-based results suggest a 
methodology for "engineering" specific two-dimensional 
patterns through a judicious choice of forcing function. 
The idea is to exploit the results in Table U in con- 
structing a multi-frequency forcing function such that 
enhancing (and/or suppressing) resonances occur at care- 
fully chosen angles. We will apply this methodology 
to demonstrate how a superlattice pattern of the SL-I 
type observed in [3| rnay be stabilized. Stabilization of 
this superlattice patterns can be related to the damped 
Q. = n — m "difference frequency" mode in Table ^ 
A demonstration is provided in llTl , and further ex- 
plorations arc performed in 19|, l2Clj| . The method we 
outline below, however, results in a dramatically more 
pronounced stabilization than was obtained in previous 
work. In particular, it can lead to stable superlattice pat- 
terns at onset of the primary instability of the flat fluid 
surface. 

Step 1: Use geometry to determine the angles for 
the desired enhancing (or suppressing) effects. For the 
SL-I pattern, the twelve dominant waves making up 
this pattern have wave vectors that lie at the vertices 
of two hexagons, one rotated by an angle 0h < 30° 
with respect to the other; see Figure ^| The stabil- 
ity of this pattern may be studied within the frame- 




FIG. 10: Schematic of the Fourier wave vectors corresponding 
to the 12 dominant waves which comprise an SL-I superlattice 
pattern. The vectors point to the vertices of two hexagons, 
one rotated by an angle Sh < 30° with respect to the other. 



work of a twelve-dimensional bifurcation problem which 
describes the competition of these superlattice patterns 
with stripes, rhombic patterns, and hexagons. This ap- 
proach is developed in [IGL l36L l37l | (the full bifurcation 
equations may be found there). A key result is that the 
stability of the superlattice pattern associated with Oh de- 
pends on coefficients in the bifurcation equations which 
we call (64, 65, fog) J where 



64 - h{eh)/\al 

65 - 6(60° - 0,O/|a| 

66 = 6(6O° + 0,O/|a| 



(51a) 
(51b) 
(51c) 



with a and h{6) appearing in ^. In particular, the su- 
perlattice pattern is favored when (I64I, I65I, |66|) are all 
small. Since h{0) may be made small in magnitude with 
"enhancing" resonances that cause spikes in 6(0), geom- 
etry dictates that we should arrange for such resonances 
to occur at one or more of the angles 6*, 60° — 6*, 60° + 9. 
For a more detailed discussion, see |l7| . 

Step 2: Use the dispersion relation and appropriate 
resonance conditions from Table Q to find a good set 
of forcing frequencies which satisfy the geometrical con- 
straints from Step 1. For our SL-I example, since we want 
to construct enhancing resonances, we turn our attention 
to the il = n — m "difference frequency" mode. We begin 
with three-frequency {m,n,p) forcing, aiming to make 
two of (1 64 1, 1 65 1, 1 66 1) small using the two difference fre- 
quency modes Vt ~ n ~ m and f2 = p — m. We choose to 
stabilize a superlattice pattern having 6h — 20.3° (this is 
a different SL-I pattern than the one observed in Q , but 
it is in the same family of patterns; see jlfll^). The two 
wave numbers corresponding to the difference frequency 
modes satisfy the resonance conditions 



n'^ikp^m) = [p-mf. 



(52a) 
(52b) 



With the optimal wave numbers for these damped modes 
dictated by geometry, the aim is to find a set of forc- 
ing frequencies {m,n,p) such that kn-m and kp-rn of 
(|52|l are as close to the optimal wave numbers as pos- 
sible. In practice we also vary Fq so as to arrange for 
frequencies {m,n,p) that are not too large - this is not 
strictly necessary but it eases our numerical computa- 
tions to use smaller sets of integers. In this case we ob- 
tain reasonable agreement by using {m, n,p) = (8, 10, 11) 
and Fo = 5.24. The wave numbers predicted by (|5^ are 
{kn-m, kp-m) — (0.351, 0.682) and the corresponding res- 
onance angles of © are {9n-m, Op-m) - (159.8°, 140.1°). 
These will cause spikes in b{d) at appropriately 20.2° and 
39.9°, respectively; note that the latter angle is close to 
60° -20.3° (c/. Eq.EIlD). 

We compute the coupling coefficient from H36|l with 
damping 7 = 0.1, forcing amplitude ratios |/n|/|/m| — 
1-54, |/p|/|/,n| = 1.85, and forcing phases (0„i, </)„, 0p) = 
(0°,0°,0°). The forcing ratios were chosen to make 61.0s, 
which is proportional to |/„P in one case and |/pp in 
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the other, as large as possible, while at the same time 
avoiding the critical values {i.e., the modes oscillating 
at n/2 and p/2 must remain damped). The coefficients 
(64,65,66) are represented in Figure [TTk . where we plot 
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FIG. 11; (a) Coupling coefficient for computing superlattice-I 
pattern stability. We use three-frequency forcing with 
{m,n,p) = (8,10,11). The two small "bumps" at 
{O„-m,0p-m) ~ (20.3°, 39.9°) are due to resonance with the 
modes oscillating with the difference frequencies Q — n ~ m 
and n = p — m. No superlattice patterns are stable, (b) 
Like (a), but with additional forcing frequency components 
{q, r) = (4, 6) which parametrically force the difference fre- 
quency modes. The result from (a) is duplicated as a dotted 
line for comparison. The two bumps become two very large 
spikes, and the superlattice pattern with angle Oh — 20.3° 
is stabilized. The fluid parameters used are 7 = 0.1 and 
Fo = 5.24. The forcing amplitude ratios and phases used are 
given in the text. The region around 60° corresponds to a 
hexagonal interaction not captured by our calculation, and 
thus has been removed. 



too weak to stabilize a pattern at the chosen angle, and 
so more work must be done. 

Step 3: Use the results in Table Q| to further en- 
hance/suppress the nonlinear interactions. In this case 
we add the forcing components (g, r) — (4, 6) in or- 
der to parametrically force the damped $7 = 2 and 
ri = 3 difference frequency modes and obtain larger |6rcs|- 
In order to favor our chosen SL-I pattern, we choose 
\fn\/\frn\ and |/p|/|/m| as before, and take \fq\/\fm\ = 
0.184 and |/r|/|/m| = 0.505. These ratios are close to 
(but below) their critical values when \fm\ = |/m|'^"*- 
We have chosen the phases to be {(()„%, 4'n,4'pi4'q^4'r) = 
(0°,0°,0°,-7°,-10°). Though the arguments of Section 
m suggest that we should make 6ros as large and positive 
as possible to favor the pattern, we are working with a 
cubic truncation of the bifurcation equations and so we 
actually want 6ics such that |6| is very small (as previously 
stated). We might have adjusted the forcing amplitude 
ratios to achieve this situation, but instead, we find it 
more convenient to vary the forcing phases away from 
the optimal values predicted by Table Q] 

The coupling coefficient appears as the solid line in 
Figure lllb . It nearly duplicates the result from the 
three-frequency case (which is included as a dotted 
line for comparison) but the two small bumps have 
become large spikes. We find that at 6* = 20.3°, 
(64,65,65) = (0.02230,-0.01887,-0.0045). To study 
the stability of the superlattice states, we perform a 
bifurcation analysis using the overall forcing strength 

hot = VI/mP + IM' + l/pl' + l/J' + IM' as the bi- 
furcation parameter. A branch of superlattice patterns 
with 9h ~ 20.3° bifurcates transcritically from the trivial 
state, and the subcritical branch then turns around in a 
saddle-node bifurcation at a particular value ftot = fsN- 
At a slightly greater forcing strength fsL > fsN (still in 
the subcritical regime), the superlattice pattern is stabi- 
lized, and remains stable for ftot > fsL (at least within 
the realm of validity of the weakly nonlinear description 
provided by the bifurcation equations). 

The methodology here is more successful than our pre- 
vious attempts at stabilizing superlattice patterns. Our 
work in [l3| created a spike at only one angle (as opposed 
to two, as here) and that in 'l9| did not parametrically 
force the damped mode and did not make an appropri- 
ate choice of phases in order to optimize the size of the 
spike. By combining multiple resonances with appropri- 
ately chosen phases, we have used TableQlto obtain dra- 
matically increased stabilization of the desired pattern. 



C. A conjecture on quasipatterns 



b{6)/\a\ as a dotted line. As expected, there are two 
bumps due to the two difference frequency resonances, 
though they are quite small (the large dip around 9 = 0° 
is due to resonance with the H = m mode). In fact, 
though the observed resonances at 20.3° and 39.9° are in 
excellent agreement with the prediction, the effect is far 



The superlattice pattern discussed above belongs to 
one intriguing class of complex patterns; another such 
class is that of quasipatterns. Quasipatterns are the con- 
tinuum analogues of quasicrystals. Unlike the superlat- 
tice patterns, they are not spatially periodic. However, 
their Fourier spectra possess discrete rotational symme- 
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try. Quasipatterns have been observed in a number of 
Faraday wave experiments, including [gIItLIssII. 

A common approach to certain types of quasipatterns 
has been to describe them using amphtude equations for 
the evolution of a number of critical modes equally spaced 
around a critical circle in Fourier space; see, for example, 
|2lL I39I lioj . Recent work in lijj elucidates the technical 
problems with this approach. The issue is that through 
nonlinear interactions, the critical modes generate other 
modes which come arbitrarily close to the critical cir- 
cle, and a center manifold reduction to a finite dimen- 
sional bifurcation problem is not possible. The usual 
amplitude equation description is thus without a rigorous 
mathematical foundation. Nonetheless, our basic phys- 
ical ideas should still apply to quasipatterns. We may 
tune our forcing function to drive energy into modes cor- 
responding to different resonant angles and thus favor the 
corresponding patterns. 

For example, here we suggest a forcing function which 
may favor a 14-fold quasipattern, which, to date, has not 
been observed in Faraday wave experiments. We use the 
methodology outlined in the previous example. We wish 
to arrange for b^-cB > at the angles 6j — j(180°)/14, 
j — 1 ... 3. This actually accounts for twelve of the an- 
gles in the pattern, since by symmetry, b{9) = &(180° — 
6) ~ 6(180° -|- 9). We choose the seven- frequency forc- 
ing function (m, n, p, 9, r, s, i) = (12, 17, 20, 27, 10, 16, 30) 
with the capillarity parameter Fq — 28.8 in (|36|l . The 
ri — n — m = 5, ^=p — m = 8, and fl = q — m = 15 
difference frequency modes are parametrically forced by 
the {r,s,t) = (10,16,30) components. We take forc- 
ing frequency ratios |/™l/|/m| = 1-2, |/p|/|/m| = 1-6, 
\U/\U = 2.8, \fr\/\fm\ = 0.62, \fs\/\f,n\ = 1-2 and 
\ft\/\fm\ = 2.2 and compute b{6) from (|36|l . The results 
are shown in Figure [T^ The dotted line corresponds to 
the naive choice of zero for all of the forcing phases. The 
solid line corresponds to the optimized case prescribed 
by Tabled namely (j)r = 4>s = 4>t = 90°. In both cases, 
the three difference frequency modes cause spikes in b{9) 
at the required angles. 



VI. CONCLUSIONS 

In this paper, we have used methods from equivariant 
bifurcation theory to study resonant triad interactions 
in Faraday waves. We have shown how the spatial and 
weakly broken temporal symmetries (or alternatively, pa- 
rameter symmetries) may be used to determine which 
spatiotemporally resonant damped modes play the most 
important roles in pattern selection. The symmetry- 
based analysis not only identifies the modes, but tells 
us how the strength of the triad interactions depend on 
the frequencies, amplitudes, and relative phases of the 
various components in an arbitrary multi-frequency pe- 
riodic forcing function. In many cases we know whether 
the interaction has an enhancing or suppressing effect on 
associated patterns. The study in this paper constitutes 
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FIG. 12: Coupling coefficient b{9) in @ as com- 
puted from l^fcil l for seven-frequency {m,n,p,q,r,s,t) = 
(12,17,20,27,10,16,30) forcing. The forcing phases are all 
for the dotted line, while the solid line corresponds to the 
optimal choice (f>r,s,t = 90°. The three damped modes with 
frequencies n — m,p — m, and q — m are parametrically forced 
by the (r, s, t) components. These three difference frequen- 
cies lead to spikes in b{0) at angles which may help stabilize 
a 14-fold quasipattern; the desired locations of these spikes 
(as determined by geometry) are indicated by dashed vertical 
lines. The small spike around 42° is due to another difference 
frequency resonance not of interest here. As in Figure ITTI the 
region around 60° has been removed. 



a somewhat unusual situation (we know of only a few oth- 
ers, such as I42J) because significant information about 
the bifurcation coefficients, namely their scaling with re- 
spect to the physical parameters and in some cases their 
sign, can be obtained without resorting to calculations 
using the governing equations. This is possible because 
of the structure imposed by the parameter symmetries of 
the problem. 

We have applied our results to impulsively forced and 
multi-frequency forced Faraday waves in several exam- 
ples, emphasizing how the resonant interactions can be 
controlled by choosing judiciously the parameters in the 
forcing function /(t). An appropriate choice allows one 
to stabilize complex patterns such as the superlattice-I 
pattern examined in Section Techniques based on Ta- 
ble^may be useful to experimentalists wishing to observe 
specific patterns in the laboratory. 

The results in this paper tie together many of the ideas 
explored in [VA Il9l 12(1 12^ and provide an exhaustive 
description of the important resonant triad interactions 
for Faraday waves (with sufficiently weak damping) . Re- 
cent experiments used multi-frequency forcing of Fara- 
day waves in order to control the transition between dif- 
ferent nonlinear states and to suppress spatiotemporal 
disorder 43). In particular, the authors of |43| apply a 
perturbing third frequency to two-frequency forced pat- 
terns near a codimension-two point and interpret their 
results in terms of the temporal parities of the dominant 
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forcing frequency and the perturbing frequency. Our re- 
sults in Table U suggest that the frequencies themselves 
(not just the parity) and the forcing phases are impor- 
tant, thus providing an alternative approach for control- 
ling patterns. 

It will be interesting to extend our work to other sys- 
tems. For example, in vertically vibrated convection, 
Boussinesq symmetry prohibits three-wave interactions 
|44| . Four- wave interactions are the important nonlin- 
ear interactions, and are the building blocks of complex 
square superlattice patterns observed in [4^ [ig • Apply- 
ing techniques similar to those developed here might yield 
insight into this pattern selection mechanism as well. 
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